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Abstract 

Within a semiclassical description of above-threshold ionization (ATI) we identify the interplay 
between intracycle and intercycle interferences. The former is imprinted as a modulation envelope 
on the discrete multiphoton peaks formed by the latter. This allows to unravel the complex 
interference pattern observed for the full solution of the time-dependent Schrodinger equation 
(TDSE) in terms of diffraction at a grating in the time domain. These modulations can be clearly 
seen in the dependence of the ATI spectra on the laser wavelength. Shifts in energy modulation 
result from the effect of the long Coulomb tail of the atomic potential. 

PACS numbers: 32.80.Rm,32.80.Fb,03.65.Sq 
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Tunneling ionization is a highly nonlinear quantum-mechanical phenomenon induced by 
intense laser pulses (> 10 14 W/cm 2 ). Electrons are emitted by tunneling through the poten- 
tial barrier formed by the combination of the atomic potential and the external strong field. 
Tunneling has recently attracted increasing interest as a probe of the atomic and molec- 
ular structure 

m. 

Tunneling occurs within each optical cycle predominantly around 
the maxima of the absolute value of the electric field. The interference of the successive 
bursts of ejected electrons reaching the same final momentum gives rise to features in pho- 
toelectron energy and momentum distribution which are markedly different from typical 
above-threshold ionization (ATI) spectra by multi-cycle lasers. This temporal double-slit 
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interference has recently been studied both experimentally [l|, |4| and theoretically [5|. On 
the other hand, the ATI peaks separated by a photon energy can be themselves viewed 
as an interference pattern formed by electron bursts repeated each optical cycle. Details 
of the interplay between these intra- and intercycle interferences have not yet been clearly 
identified and analyzed, to the best of our knowledge. 

In this Rapid Communication, we study the influence of different interference processes 
on ATI spectra generated by multi-cycle laser pulses. We clarify the underlying mechanism 
within a simple one-dimensional model employing classical trajectories. Within the frame- 
work of the strong-field approximation (SFA) [6( the qualitative features, the modulation 
of the ATI peaks akin to the modulation of Bragg peaks by the structure factor in crys- 
tal diffraction, can be unambiguously identified in the ATI spectrum determined from the 
full solution of the three-dimensional time-dependent Schrodinger equation (TDSE). The 
multi-cycle laser pulse thus acts as a diffractive grating in the time domain. Quantitative 
deviations between the SFA predictions and the full TDSE can be traced to the Coulomb 
tail of the atomic potential affecting this modulation. The latter opens up the opportu- 
nity to observe effects of the atomic potential in easy-to-obtain photoelectron spectra after 
ionization by multi-cycle laser pulses. 

Our simple semiclassical model of photoelectron spectra is based on the ID "simple man's 
model (SMM)" . Let us consider an atom interacting with a linearly polarized laser 

pulse. The laser field F(t) is chosen to be of the form F(t) = f(t) smut, with an envelope 
function f(t) corresponding to an iV-cycle flat-top pulse with a field strength of F Q and with 
m-cycle linear ramp-on and -off (we set m = \ in the following). The classical electron 
trajectories i (i = 1, 2N) for a final momentum k are characterized by their release times 
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tr which satisfy, 

k = -A{tf\ (1) 

where A(t) = — f_ . dt'F(t') denotes the vector potential divided by the speed of light, 
with A(t) = (Fq/u) cos out for the flat-top segment of the pulse. In this study we focus on 
direct photoelectrons (without rescattering) with energies E < 2U P (U p = Fq/Auj 2 ) which 
dominate the total ionization probability. It should be noted that within the flat-top part of 
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is a constant 



the pulse and for a given value of k, the field strength upon ionization 
independent of the release time t z r (i = 1,...,2N). Thus, assuming that the ground-state 
depletion is negligible, the ionization rate r(^) is identical for all the ionization bursts 
(or trajectories) to first approximation and only a function of T(k). Consequently, the 
momentum distribution P(k) can be written as 



P(k) = T(k) 
where S denotes the Volkov action 
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with I p being the ionization potential. The key to the analysis of intracycle and intercycle 
interferences is that the sum over interfering trajectories (Eq. [2]) can now be decomposed 
into those associated with two release times within the same cycle and those associated with 
release times in different cycles (see Fig. 1) 

2N N 2 
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Accordingly, Eq. ([2]) becomes 

P(k) = r(Jfe) cos 2 



AS 



N 
3=1 



T(k)F(k)B(k) 



(5) 



where Sj = [S^) +S(t?'* } )]/2 is the average action within each cycle and AS" = S(tr 
S(tr' 2 ^) is the difference. Note that AS is independent of j. Eq. ((Sj) is structurally equivalent 
to the intensity for crystal diffraction: the factor F(k) = cos 2 (AS/ 2) represents the form 
factor (on structure) accounting for interference modulations due to the internal structure 



within the unit cell while the second factor B(k) 



gives rise to Bragg peaks 
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FIG. 1: (color online) Electric field F(t) (left axis) and vector potential A(t) (right axis) of a 
flat-top pulse (N = 4, m = |). The electron emission times for a given final momentum k are 
marked by circles (tf? ) and triangles (tr )■ Each circle-triangle pair determines the structure 
factor F(k) and leads to intracycle interference while the periodic train of such pairs gives rise to 
intercycle interference. Each optical cycle can be viewed as "unit cell" of the time lattice. 

due to the periodicity of the crystals. Alternatively, Eq. ([1]) may be viewed as a diffraction 
grating in the time domain consisting of N slits and with F(k) being the diffraction factor 
for each slit. Interferences between different slits correspond to intercycle interferences while 
F(k) represent intracycle interferences. 

We show in the following that the interplay between B(k) and F(k) controls the direct 
ATI spectrum (Fig. 2). Note that the peak of the intracyle structure factor F(k) are not 
equispaced in energy (shown in Fig. 2 (a) for hydrogen with Ip = 0.5 a. u.). The separation 
of consecutive peaks is larger at intermediate energies than near the classical boundaries 
E — and E = 2U P (U p = 0.5 a.u. in the present case). It can be analytically shown that 
the separation between adjacent peaks Q(k) as a function of the final momentum is given 
by 

m = wr=g_ 

1 — k 2 + 7 2 /2 — Ky 1 — k 2 arccos k 
where 7 = ^/2I p uj / Fq is the Keldysh parameter, and k = uk/F Q is the scaled dimensionless 



momentum. Q(k) reaches a maximum at n m = y 1 + ^ — |a/2 + 7 2 , the position of which 
is indicated by a vertical arrow in Fie. [2] (a). Such intracycle interference patterns have 
been recently experimentally observed for ultrashort near-single cycle IR pulses |4| and 
theoretically analyzed 5j. The intercycle interferences pattern (Fig. 2b) is, by contrast, 
equispaced and corresponds to the well-known ATI pattern with peak spacing corresponding 
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FIG. 2: Semiclassical intracycle and intercycle interference pattern (color online), (a) Intracycle 
interference pattern given by the structure factor F{k). (b) Intercycle interferences with "Bragg" 
peaks given by B(k) in Eq. ([7]). (c) Total semiclassical interference (Eq. ([5])) for N = 2. The 
laser wavelength and intensity are 1000 nm and 1.6 x 10 14 W/cm 2 , respectively. To stress the 
interferences, we set T(k) = 1 and normalize the respective maxima of F(k) and B{k) to unity, see 
Eq.©. 

to one-photon transitions (huj). For the case of N = 2 in Fig. 2(b), B(k) is given by 

fs(tP)-SU 1) + 2n/u J )\ 
B{k) = 4 cos 2 £ L . (7) 

In multi-cycle photoelectron spectra with N > 2, both intracycle and intercycle interfer- 
ences are simultaneously present. In Fig. [2] (c), the resulting energy distribution calculated 
from the time grating (Eq. (jS])) is shown for N = 2. The multiphoton peaks are modu- 
lated by the intracycle interference structure factor. The ATI peaks become narrower as the 
number iV of optical cycles increases approaching, in this way, 5-peaks for infinitely long 




FIG. 3: (color online) Semiclassical total ionization probability as a function of the laser wavelength 
A and the ATI order n (see text). Shown is the spectrogram for both a one-cycle pulse (green stripes) 
and and a four-cycle pulse (red islands). Fq = 0.0675 a.u. (1.6 x 10 14 W/cm 2 ). 

pulses. On the contrary, the intracycle modulation is independent of the number of cycles. 
Thus, the sub-cycle interference, previously studied with near-single-cycle ultrashort pulses, 
is embedded and visible in ATI spectra for multi-cycle pulses, a feature apparently up to now 
not fully recognized. This effect becomes more transparent when we study the parametric 
variation of the photoelectron spectrum. In Fig. [3] we show the spectrum expressed in units 
of the (new) photon number n = (E + I p + U p )/u, as a function of laser wavelength A. 
For N = A cycles the horizontal stripes in this two-dimensional interferogram peaking near 



6 



integer values of n represent ATI peaks due to intercycle interference. Superimposed are 
tilted stripes controlled by the intracycle interference, which are also visible for the N = 1 
cycle pulse. We note that a very similar two-dimensional interferogram has recently been 



found in atom-surface diffraction 
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In view of the fact that the present intracycle and intercycle interference structures are 
derived from a simple one- dimensional SFA model neglecting realistic features of the atomic 



potential, we test its predictions against a full numerical 3D TDSE solution for hydrogen. 
The resulting interferogram for photoelectrons (Fig. 4) agrees qualitatively remarkably well 
with the one- dimensional SFA model (Fig. 3). The intracycle modulation is best seen when 
the angular acceptance of the photoelectron spectrum is restricted to a cone of small angles 
around the polarization axis (9 = 10° in Fig. 4a) while it is somewhat blurred but still visible 
in the total spectrum (Fig. 4b). We found similar interference patterns for rare gases such as 
argon, though the details depend on the atomic potential. Projecting the two-dimensional 
distribution onto the A axis (Fig. 4c) results in a regular modulation pattern on a fine A 
scale which can be traced back to the combined effect of inter- and intracycle interferences. 
This oscillations closely resemble those previously observed for the A dependence of HHG 

Experimentally, more easily accessible is the variation of the photoelectron spectra as 
a function of the energy at fixed wavelength, i.e., projection of the 2D interferogram onto 
the A axis (Fig. 5). The photoelectron spectrum clearly reflects the modulation of the 
regularly spaced ATI peaks (intercycle interference) by intracycle interference, in complete 
agreement with the prediction of a time grating (Eq. (jSJ). It is now highly instructive to 
probe this interference pattern for its dependence on the Coulomb potential neglected in 
the SFA derivation of the time grating. Within the full TDSE we can explore Coulomb-tail 
effects by using a screened Coulomb potential, 

V(r)={"' ( „ <r<rc) , (8) 

( r ( r > Tc > 

where the parameter characterizes the width of the transition from the full to the screened 
Coulomb potential. We use = 10 a.u., and r c is varied between 10 and 70 a.u. For these 
parameter values, the ionization potential and the first excitation energy remain unchanged 
to an accuracy ranging between 10 -3 and 10~ 9 . Consequently, the position of the ATI peaks 
remains unchanged as a function of r c . However, the modulation envelope is systematically 
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FIG. 4: (color online) Two-dimensional interferogram as a function of laser wavelength A and 
ATI order n (see text) calculated by using the TDSE for hydrogen for a four-cycle pulse with 
F Q = 0.0675 a.u. (1.6 x 10 14 W/cm 2 ). (a) Emission into a cone of 10° around the polarization axis 
and (b) for all angles, (c) Energy-integrated total ionization yield as a function of A. 

and almost rigidly shifted towards lower energies under the influence of the long-range 
Coulomb tail. The Coulomb phase does not affect the Bragg peaks but changes the structure 
factor F(k) of the time grating in Eq. (jSJ). Thus, the Coulomb potential is responsible for 
a shift in the positions of the intracycle interference stripes, visible in Figs. |3] and H](a). 
It should be emphasized that the shift is equally significant even for higher photoelectron 
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FIG. 5: (color online) Comparison of the photoelectron spectra within the 10° cone, calculated with 
the full Coulomb potential (marked as "full") with the truncated Coulomb potentials for varying 
values of r c as indicated. The laser wavelength is 1000 nm. Other laser parameters are the same 
as in Fig. SJ Thin solid lines are approximate peak envelopes. 

energies, for which one might expect the SFA to be better justified. We note parenthetically 
that such shifts can be partially accounted for by including the Coulomb- Volkov phase into 
the semiclassical theory Q, [l5j. 

In summary, we have identified the interplay between the intra- and intercycle interfer- 
ences of electron trajectories in ATI spectra by intense multicycle laser pulses. The former, 
which carries information on attosecond subcycle dynamics of the electron cloud [lj], is im- 
printed as a pulse-length independent modulation of multiphoton peaks formed by the latter. 
This modulation is even more clearly visible in the wavelength dependence of ATI spectra, 
which could be measured with state-of-the-art tunable sources based, e.g., on optical para- 
metric chirped pulse amplification. While the general features of the interplay between the 
intracycle and intercycle interferences is well explained by the simple semiclassical theory, 
the intracycle modulation envelope is shifted by the Coulomb tail which emphasizes the 
effect from the atomic potential. 
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